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1.  INTRODUCTION  AND  SUMMARY 

In  a  previous  report  [l],  the  use  of  HAD  (mean  absolute  deviation) 
in  the  analysis  of  demand  variability  as  proposed  by  R.  G.  Brown  [3] 
and  presently  used  by  NavSup  [4]  was  discussed.  In  particular,  the  use 
of  exponentially  smoothed  estimates  of  HAD  was  criticized  on  theoretical 
grounds  and  a  case  was  made  for  further  study  into  the  probability 
distributions  of  the  various  estimators  presently  in  t.se.  Of  special 
Importance,  the  use  of  MAD  in  setting  safety  levels  was  mentioned  as  a 
problem  area  in  need  of  further  investigation. 

The  present  report  concerns  itself  with  the  extension  of  research 
investigations  in  the  areas  just  mentioned.  Several  attempts  to  derive 
the  necessary  probability  distribution  theory  for  MAD  estimators 
resulted  in  disappointing  futility  for  the  most  part.  This  was 
anticipated  for,  as  previously  mentioned  in  [1],  absolute  deviation  as 
a  measure  of  variability  has  been  abandoned  for  years  on  theoretical 
grounds  by  Statisticians. 

Two  minor  results  that  were  obtained  concern  the  ratio  of  MAD,  A, 
to  standard  deviation,  o.  The  complete  result  of  this  ratio  has  been 
determined  for  the  Poisson  case.  As  alluded  to  in  [1],  the  universal 
use  of  fl  *  0.8  for  this  ratio  is  especially  bad  for  low  demand 
Poisson  item.  The  same  general  criticism  can  be  made  in  the  negative 
binomial  case.  However,  it  was  found  that  for  the  choice  of  parameters 
presently  being  used  by  NavSup,  the  approximation  of  this  ratio  by 
0.8,  is  a  good  one. 
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Failing  to  obtain  the  required  distribution  theory  to  study  the 
behavior  of  NAD  estimates,  simulation  was  used  to  compare  smoothed 
estimates  of  a  to  the  usual  maximum  likelihood  estimates  mostly  tor 
the  case  of  normal  demand.  Parameter  choices  that  reflect  high  demand 
items  were  selected  in  order  to  give  HAD  the  best  possible  advantage  in 
the  comparison  with  alternate  methods  of  estimation.  Such  parameter 
choices  were  made  after  consultation  with  NavSup  personnel  in  order  to 
duplicate  situations  that  actually  exist;  moreover,  large  amounts  of  real 
data  are  available  for  such  cases.  In  every  single  case  studied,  the 
sample  variance  of  the  smoothed  estimate  of  NAD  and  consequent  estimate 
of  a  was  roughly  twice  that  of  the  maximum  likelihood  estimates  of  o. 
The  same  was  true  in  a  slow  mover  case  or  two  that  was  examined.  This 
confirms  in  a  general  way  the  observation  made  in  a  similar  study  carried 
out  by  Asher  and  Vallace  [2]  in  which  they  found  NAD  to  be  significantly 
less  efficient  (about  20%)  compared  to  the  classical  minimum  variance 
estimators  under  the  Gauss-Narkov  assumptions. 

Mot  content  merely  to  confirm  this  observation,  however,  the  estimates 
found  were  used  further  to  see  what  the  actual  effects  on  reorder  levels 
were.  This  was  approached  in  several  ways.  First,  NAD  and  maximum 
likelihood  were  compared  by  computing  the  percentage  of  time  the  known 
theoretical  reorder  level  was  overestimated,  thereby  resulting  in  too 
much  on  hand  stock.  Secondly,  the  results  were  examined  to  see  what 
percentage  of  the  observations  fell  with In  k  units  of  the  theoretical 
order  level  for  k  ■  1,2,"*  and  various  fixed  values  of  risk.  Finally, 
the  two  methods  were  compared  by  setting  reorder  levels  and  then  computing 


the  actual  risk  attained  by  those  levels  where,  for  each  cast  again,  the 
theoretical  or  true  risk  is  known.  As  seen  in  Section  4,  exponentially 
smoothed  mean  and  MAD  came  off  se'cnd  best  compared  with  maximum  likeli¬ 
hood  estimates  of  p  and  c  in  every  single  case  examined. 

Recognizing  the  limitations  of  simulation,  no  svaaping  claim  is  here 
made  for  proof  that  exponentially  smoothed  MAD  is  an  inefficient  method 
of  accounting  for  demand  variability.  At  the  same  time,  one  cannot  ignore 
the  fact  that  present  methods  were  uniformly  inferior  in  the  situations 
examined. 

In  Section  5  recommendations  for  further  study  are  made.  Among  these 
is  the  suggestion  that  real  data  be  used  from  the  histories  available  at 
FMSO  for  comparing  the  results  in  retrospect  with  what  would  have  been  the 
case  had  maximum  likelihood  procedures  been  used.  The  writer  wishes  to 

H 

thank  Lt.  Ozden  Orneck  for  his  invaluable  assistance  in  constructing  and 
running  the  computer  programs  as  wall  as  assistance  in  mathematical 
derivations.  Acknowledgement  should  also  be  given  to  Cdr.  Jack  E.  White  of 
FMSO  for  his  unfailing  cooperation  in  defining  the  problem  areas  and  supply¬ 
ing  parameter  values  that  are  realistic  in  terms  of  NavSup  use.  Credit 
should  also  be  given  to  Mi.  Jamas  W.  Prichard  (SUP  04E)  for  his 
continued  endorsement  and  interest  in  this  research  area. 

2.  RATIO  OF  A  TO  a 

It  will  be  helpful  to  review  some  parametric  definitions  and 
establish  a  notation  to  be  used  here  and  in  ensuing  sections.  Let  X 
be  a  random  variable  with  mean  p  and  variance  cr*  (standard  deviation 
u) .  The  mean  absolute  deviation,,  MAD  for  shot*  will  be  denoted  A 


and  is  defined  as 


(2.1)  A  *  E(|*n|) 

As  remarked  in  [l],  it  would  be  more,  rational  to  define  MAD  as  E(|x~m|) 
where  a  is  the  median  of  X  bar.  lit.  is  (2.1)  that  is  used  by  Brown  [3] 
and  NavSup  [4].  Consequently  the  same  definition  is  adopted  in  this 


report. 

Now  whe.,  X  is  normally  distributed It  is  well  known  that  the 
ratio  of  A  to  a  is  given  by  cr?  roughly  0.8.  It  is  somewhat 
surprising  that  the  ratio  is  approximately  the  same  for  certain  other 
families  of  probability  distributions  such  as  the  Exponential,,  Uniform 
and  Triangular.  However,,  this  result  is  not  universally  true  and  is 
particularly  a  poor  approximation  for  the  P®iss@n  family,,  a  model  often 
used  for  the  so-called  slow  mover  type  of  inventory  item.  This  fact 
was  demonstrated  earlier  in  [l]  by  examining  the  ratio  for  selected 
values  of  the  Poisson  parameter  „  X.  That  ratio  has  now  been  determined 
for  all  values  of  X  and  It  may  be  instructive  to  see  the  behavior  of 
this  ratio  in  a  complete  sense. 

By  [z],  we  shall  mean  the  greatest  integer  In  r,  that  Is,  the 
largest  integer  n  such  that  n  *  s  (and  hence  %  <  n  +  i) .  We 


denote  the  Poisson  (X)  mass  function  by  p  so  that , 


~k  X 


X  -  0..1  IB;,'" 


Then  we  have,  recalling  that  p  *  -  X  for  this  case., 

•  [X]  • 

A  -  B(|x-v!)  -  ^  |x-X|p(x)  »  ^  (X'x)p(x)  +  )  (x-X)p(x) 

x-0  X"0 

z 

But,  letting  F(z)  ■  Ip(x)  define  the  Poisson  distiibution  function, 

x"Q 


,  4  ... 


[X] 

I 


[_xj 

(X-x)p(x)  -  X  P([X])  -  ^  xp(*)  "  X  F([X] 


*-ra+i 


xp(x)  -  X 


Y  (x-X)p(x)  -  V  xp(x)  -  X  +  X  F([X]) 
x-[XJ+l  x-[T]+l 

Adding  these  results,  a 

(2.2)  A  -  2  Xp([X])  -  2  X  +  2  V  xp(x)  . 

x"[XJ+l 


But,  for  x  1  1, 

-X  ix  -x  ,x-l 

Xp(x)  -  xe  ■  Xe  X  ■■■-;■  -  Xp(x-l)  . 
x.  (x-1); 

Hence, 

_•  •  • 

y  *p(x)  *  x  y  p(x-i)  ■  x  y  p(y>  -  x(i-?([x]-i)). 

x"[X]+l  x"[Xj+l  y"[X] 

Substituting  this  result  in  (2.2), 

A  -  2  X  F([X])  -  2  X  +  2  X  -  2  X  F([X]-1) 

-  2  X(F[X]-F([X]-1)) 

-  2  X  p([X]) 


Thus, 


(2.3) 


A  ■  2  Xe 


-X 


iw 

njr 


for  any  X  >  0. 


Since  a  •  / 1  we  have, 
(2.4) 


Cx]: 


The  graph  of  “  as  a  function  of  X  appears  in  Figure  1. 


^  « 


0 


1 


2 
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Figure  1.  Graph  of  - 

The  analysis  that  validates  the  graph  is  the  following. 

If  n<X<n+l  then  [X]  ■  n  and,  ^  ■  -|i  X  e  . 

As  a  function,  say  f  ,  of  X  in  this  open  interval,  f  is 

n  n 

differentiable  with, 

f'(X)  -f.  «‘X  X“'4 

n  n* 

Hence,  f  #(X)  >  0  if  X  <  n  +  %  and  f  *(X)  <0  if  X  >  n  +  £ 
n  n 

so  that  fn  possesses  a  relative  maximum  at  X  ■  n  +  \  with  a 

maximum  value  of  fQ(n+4)a  for  n  ■  0,1,2,***. 

With  the  formulas  established  for  f  and  f  ,  it  is  easy 

n  n+1 

A 

to  verify  that  ~  is  continuous  at  each  of  the  integers  although 

not  differentiable  there.  In  any  case,  it  is  shown  in  the  appendix 

A 

that  the  sequence  of  values  of  -  at  n, 
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is  monotone  increasing.  Using  Stirlings  approximation  [7]  to  n!, 

m  nn+1  .•”<»:<  m  .'"(i+tt-t) 

12n-l 

ve  see  that, 

fl - -  <  £  (n)  <  fl 

tt  l+~~-  tt 

I2n-1 

and  lim  £  (n)  ■  fl  .  Thus  the  values  of  -  converge  in  the  limit 
n-»»  Q  ir  0 

to  the  factor  fl  .  However,  as  previously  remarked  in  [l]  and  as 
TT 

seen  vis mlly  in  Figure  1,  approximating  ^  by  Jl  is  particularly 

a  TT 

bad  when  X  <  1  which  is  precisely  when  the  Poisson  assumption  is 
of  primary  interest. 

For  intermediate  movers,  the  negative  binomial  distribution  is 
employed  as  a  model  for  demand.  Specifically,  from  [4],  such  an 
assumption  is  made  whenever  mean  demand  is  between  2  and  20  units. 
The  negative  binomial  distribution  is  a  family  of  probability 

i 

I 

mass  functions  having  two  parameters  r  and  p  where  r  is  a 
positive  integer  and  0  <  p  <  1  .  For  our  purposes  the  general 
formula  for  the  mass  function  is  defined  by, 

p(x)  -  C1*”1)  pV  ,  x  -  0,1,2,* **  (q=l-p) 

The  mean  li  and  the  variance  of  a  random  variable  having  this 
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mass  function  are  given  by, 


I.-S  ,  o2--*? 

P  ’  p2 


Proceeding  as  in  the  Poisson  case, 

»  DO  » 

i  ■  l(|x-u|)  -  ^|x-p|p(x)  -  ^(p-x)p(x)  +  V  ^(x-p)p(x) 


M  Cp] 

1'(P-x)(^SpV  -  ppr  +  Ifc-xH^-bpV 

x"0  x"l 


LW 

I  (x-PM^-^V  -  -  Ix^W  +  P I  i"?W 

x-fiTj+l  x-1  X  x*0  * 

Adding  these  results, 

M  Cp] 

(2.5)  A  -  2[ppr  +  p^p(x)  -  ^xp(x)] 

expressing  A  in  terns  of  finite  suns. 

y~rq  A 

Since  a  ■  “  ,  it  is  no  problem  to  compute  the  ratio  -  for  various 

choices  of  r  and  p  on  the  computer. 

It  is  easy  to  choose  values  of  r  and  p  for  which  the  ratio 

is  very  different  from  0.8.  However,  the  interest  would  only  be 

academic  for  values  of  r  and  p  such  that  (i  {:  [2,20].  For  that 

range  of  values,  a  total  of  1,000  cases  were  considered  as  follows. 

Let  €  [2,20]  be  selected  and  for  each  such  r  was  allowed 

to  vary  over  the  set  {2,4,6,*** ,20} .  For  each  such  choice,  p  is 

r 

then  determined  from  the  relation  p  ■  —  and  q  is  specified  by 

r  -H* 


q  ■  1  *  p.  Then  A  was  computed  using  (2.5)  and  finally  “  was 
evaluated.  In  each  such  "group"  of  ten  parameter  pairs,  it  was 
observed  that  ~  was  invariably  monotone  Increasing  although  there 
was  no  monoticity  from  group  to  group.  More  significantly  however, 
it  was  observed  that  -  varied  always  between  a  minimum  value  of 
0.750  and  a  maximum  value  of  0.805.  For  example,  if  p  ■  0.893 
and  r  ■  20  whence  ^  ■  2.4,  we  find  that  £  B  0.805  although  for 

p  ■  0.455  and  r  *  2,  resulting  also  in  ■  2.4,  ~  ■  0.770.  On 

the  other  hand,  for  p  ■  0.174  and  r  ■  4,  yielding  ^  *  19, 

j  ■  0.781  so  that  low  values  of  g  are  not  always  associated  with 

low  values  of  i*. 

In  summary,  then,  we  may  say  that  for  the  negative  binomial 
case  and  the  range  of  parameter  values  presently  used  by  NavSup, 
the  approximation  /2  for  ^  appears  to  be  a  safe  one. 

3*  ESTIMATING  VARIABILITY  WITH  MAD 

Having  established  the  value  of  ^  for  cases  of  particular 
interest  in  the  inventory  models,  it  follows  that  o  ■  k  A  for 
some  constant  k.  Consequently,  if  A  is  estimated  directly  from 

aw 

the  data  to  obtain  an  estimate  A  ,  say,  then  o  can  be  estimated 

m  aw 

by  applying  the  formula  a  ■  k  A  .  In  this  way,  one  obtains  a 
measure  of  the  variability  of  demand  through  estimates  of  A  . 

How  good  such  a  procedure  might  be  clearly  depends  upon  the 
technique  used  to  estimate  A  in  the  first  place  and  what  the 

M 

effect  is  on  scaling  that  estimate  by  the  factor  k  to  obtain  o. 
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There  ere  of  course  many  different  ways  to  estimate  o  . 

Following  Brown  [3],  NavSup  uses  a  formula  based  upon  exponentially 
smoothing  certain  error  forecasts  defined  as  follows*  Let  XQ,Xjy  *,Xt 
denote  independent  demands  over  discrete  units  of  time  up  to  and 
including  t.  Suppose  that  £  ^  is  a  forecast  of  the  demand  at 
time  t,  such  forecast  being  made  at  time  t-1.  Then  the  difference 
et  ■  Xj.  -  ,  represents  a  forecast  error  that  is  dependent  upon 

the  method  of  forecasting  used.  Again  following  Brown,  suppose  that 
exponential  smoothing  is  used  to  forecast  so  that,  using  smoothing 
constant  Of, 


t-2 


In  this  report  we  will  only  be  considering  the  model  where  random 

2 

demand  X  is  normal  with  constant  mean  )i  and  variance  o  .  With 
this  in  mind  we  have,  as  previously  reported  in  [l], 

S  (et)  ■  0  for  all  t 

and  the  error  terms  have  time ''dependent  variances  given  by, 


o\(t)  -  a2  +  o2  <£|£ 


2t“l 


However,  if  we  let  t  ■*  *  ,  we  obtain  a  limiting  variance  of 


(3.1) 


2  2  2 
<J  -  -f-  O  or  o 

«  2-0 


/T~a 

2 


a 

e 


Now  for  each  t,  f  ,  being  a  linear  combination  of  independent 
normal  random  variables,  is  again  normal  and  its  mean  absolute 
deviation,  Ae,  is  related  to  by  ae  “  Ae  *8  Previ°u8^y 


C3 


10 


established.  Consequently,  from  (3.1)  ve  nay  write,  asymptotically, 


(3.2) 


This  formula  relates  demand  variability  to  the  mean  absolute 
deviation  cl  error  forecasts. 

The  formula  used  in  [3]  and  [4]  to  estimate  A  and  ultimately 
o  Is  based  on  an  exponentially  smoothed  function  of  the  error 
forecasts.  If 


then  o  Is  estimated  for  demands  by 


(3.3) 


~  /ir(2-oft 

a"  2 


A 

t 


As  remarked  in  [l],  the  worth  of  this  estimation  procedure  is 
difficult  to  assess  without  some  knowledge  of  the  probability 
distribtuion  of  A^  .  Serious  efforts  notwithstanding,  that 
probability  distribution  has  not  been  derived  to  date.  The  real 
basic  difficulty  is  that  random  variables  ei*«2»***»et 
compose  At  ,  while  normally  distributed,  are  not  mutually 
independent.  Even  the  matter  of  finding  the  second  moment  of  A^ 
has  been  intractible  thus  far. 

In  an  effort  to  discover  some  indication  of  the  variability 
involved*  a  simulation  study  was  undertaken  and  resulted  in  mounting 
evidence  that  there  is  reason  to  believe  that  the  variance  of  A^ 
is  bounded  below  by  a  positive  quantity  so  that  even  in  the  limit 


•=•  11 


the  actual  estimates  o  based  on  this  quantity  will  fluctuate 
about  the  true  value  a. 

In  turn,  this  immediately  suggests  comparing  this  method  of 

A 

estimation  with  the  classical  maximum  likelihood  estimate  a 
given  by  the  formula, 

(3.4)  a  - 

t 

Now,  properties  of  o  are  well  known  (see  [5]  for  example) ,  being 
derived  from  the  so-called  Chi-distribution.  In  fact, 

K(a)  .fl  o^a  where 

and, 

V(o)  ■  pj  a2  ,  where  ■  1  -  cr2  *  — “ — 
t  t  t  2(t-I) 

From  the«e  formulas  we  see  that 

-  a 

(3.5) 

lim  V(o)  ■  0 

t400 


Hence,  if  it  is  true  that  the  variance  of  o  is  bounded  away  from 

A 

zero,  then,  even  though  unbiased  in  the  limit,  the  fact  that  o  has 
a  vanishing  variance  would  make  it  preferable  as  an  estimating  tool. 
Our  simulation  results  certainly  seem  to  concur  that  this  is  the 
case. 

For  simulation  purposes,  several  members  of  the  normal  family 
were  selected  so  as  to  be  representative  of  a  wide  class  of  fast 
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moving  Inventory  items.  Using  random  generation  samples  of  size 
200  were  drawn  from  these  populations  and  a  running  account  was 

+•  A 

kept  of  the  estimates  o  and  o  of  the  true  and  known  standard 
deviation  9.  This  was  done  by  computing  the  values  of  (3.3)  and 
(3.4)  at  each  time  period  t. 

Incidentally i  it  should  be  observed  in  this  regard,  that  the 

A 

running  computation  of  9  involves  no  more  time  nor  storage  than 

that  of  9,  a  feature  often  cited  as  one  compelling  reason  for 

*2 

using  exponential  smoothing.  It  is  true  that  a  will  have  to  be 
stored  and  updated  at  the  next  time  period  and  hence  one  extra 
operation,  that  of  taking  the  square  root,  will  be  involved.  But 
with  modern  computers ,  the  time  for  this  extra  operation  is 
negligible. 

The  actual  point  by  point  results  of  the  simulation  are  perhaps 
not  too  enlightening.  The  program  that  was  used  has  been  preserved 
and  is  available  for  further  use.  Of  more  Interest  is  a  comparison 
of  the  two  procedures  with  regard  to  bias  and  vai  iance.  These 

have  been  estimated  by  computing  the  sample  average  and  sample 

*  •* 

standard  deviation  (S.D.)  of  each  of  a  and  9  for  each  parameter 
choice.  While  many  more  cases  were  examined ,  the  results  for 
parameter  choices  furnished  by  NavSup  are  summarized  in  Table  1. 


Parameter 

Pairs 

Average 

aw 

a 

Average 

A 

a 

Bias 

« 

a 

Bias 

A 

a 

S.D. 

a 

m 

M.S.K. 

a 

m 

(50,10) 

10.27 

9.60 

0.27 

2.83 

1.33 

8.08 

1.93 

(100,25) 

25.67 

23.99 

0.67 

7.08 

3.31 

50.57 

11.98 

(100,50) 

51.33 

47.98 

1.33 

14.15 

6.63 

201.99 

48.04 

(400,40) 

41.07 

38.39 

1.07 

-1.61 

11.32 

5.30 

129.29 

30.78 

(500,50) 

51.38 

48.52 

1.38 

=1.48 

14.15 

7.40 

202.13 

56.95 

(600,60) 

61.60 

57.58 

1.60 

-2.42 

16.98 

7.95 

290.88 

69.06 

(700,70) 

64.99 

58.78 

■5.01 

-11.22 

21.40 

10.56 

483.06 

237.40 

(800,80) 

81.51 

76.15 

1.51 

-3.85 

19.18 

5.85 

370.15 

49.04 

(900,90) 

87.26 

86.21 

*2.74 

-3.79 

21.09 

5.01 

452.30 

39.46 

(1,000,109) 

101.89 

95.18 

1.89 

578.13 

76.81 

Table  1.  Sample  Characteristics 


The  table  clearly  brings  out  the  inefficiency  of  the  smoothing 
technique  compared  with  classical  methods  of  estimation.  Both  proce¬ 
dures  are  biased.  The  smoothing  procedure  tends  to  overestimate  while 
the  maximum  likelihood  procedure  more  conservatively  underestimates 
o,  a  fact  known  from  the  theory  of  course.  Bxcept  for  the  one 

A 

extreme  case  (the  pair  (700^70)  where  the  bias  for  o  is  numerically 
as  large  as  11  for  some  reason  that  is  not  consistent  with  the 
other  results,  the  bias  in  both  cases  is  of  the  same  magnitude.  There 
is  a  marked  difference  in  the  variance  estimate  of  the  two  procedures, 

M 

however.  The  variance  for  o  is  significantly  nigher  than  that  for 

A 

a  in  every  case.  Indeed,  there  was  never  a  single  normal  case 
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studied  in  which  this  was  not  the  case. 


Since  the  bias  is  roughly  the  same  for  each  procedure  this  means 
that  the  mean  squared  error  (M.S.E.)  of  each  procedure  will  largely 
be  determined  by  its  variances.  Recalling  that  M.S.E.  is  given  by 
the  variance  plus  the  square  of  the  bias  we  use  the  averages  and 
standard  deviations  of  Table  1  to  compute  estimates  of  M.S.E.  which 
are  recorded  in  the  last  two  columns.  Again  the  difference  in  the 
tv/o  procedures  are  quite  striking.  The  M.S.E.  for  o  reaches  as 

A 

high  as  12  times  that  for  0  in  one  case  and  even  in  the  pathological 

A 

case  cited  above  (the  pair  (700,70))  where  the  bias  of  a  is 

w  A 

unexplainably  high,  the  M.S.E.  of  a  is  still  twice  that  of  o. 

In  most  of  the  cases,  the  ratio  is  about  4. 


A  similar  analysis  was  undertaken  for  Poisson  demand  with  results 
that  are  not  quite  as  striking.  These  are  sunnarized  in  Table  2  for 
three  typical  cases. 


i  Parameter 

H 

Average 

0 

Average 

A 

a 

Bias 

mt 

a 

Bias 

A 

0 

1 

B 

M.S.E. 

0 

M.S.E. 

s 

j 

.01 

0.1 

.01 

.06 

O 

• 

9 

o 

• 

• 

.04 

.04 

.0032 

0.1 

0.3 

.17 

.23 

-.13 

-.07 

.13 

.0085 

1.0 

1.0 

.91 

.95 

1 

• 

O 

VO 

k 

• 

O 

<«n 

.28 

.12 

.0169 

Table  2.  Poisson  Demand 

Again  it  should  be  observed  that  the  M.S.E.  for  o  is  at  least  3 

A 

times  that  for  o  in  each  instance  examined. 
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4.  EFFECTS  ON  REORDER  LEVELS 

Not  satisfied  to  merely  summarize  the  evidence  obtained  as  in 
the  tables  of  the  preceding  section,  it  was  deemed  advisable  to  test 
the  effects  of  the  forecasting  errors  on  their  use  in  fixing  reorder 
levels.  Again  for  reporting  purposes,  the  main  concentration  was  spent 
on  the  normal  case.  If  demand  is  normal  with  mean  and  standard 
deviation  o,  then  the  (theoretical)  reorder  level  would  be  set  at 
|i  +  ko  where  k  is  chosen  to  satisfy  a  given  risk  p  defined  by 

p  ■  P(X  >  +  ka) 

Thus,  the  risk,  for  present  purposes,  is  the  probability  that  demand 
will  exceed  the  reorder  level  causing  a  stockout.  For  given  p  (or  k) 
the  parameter  k  (or  p)  can  be  determined  from  standard  normal  tables. 

Of  course  p  and  a  are  unknown  so  that  p  can  never  be 
satisfied  exactly.  If  we  were  to  use  and  (J/  for  the  reorder 
level  and  it  should  happen  that  +  ko/  >  p.  +  ka  where  k  has 
been  chosen  to  satisfy  a  given  risk  requirement  p,  then  the  true 
risk  say  p*  ■  P(X  >  +  ko*)  would  be  something  smaller.  This 

means  that  simultaneously  overestimating  n  and  o  results  in  over¬ 
stocking,  that  is,  we  could  stock  less  to  achieve  the  required  risk 
level.  There  may  also  be  penalty  costs  in  the  way  of  storage  cost6 
to  consider  for  such  a  situation.  On  the  other  hand,  if  +  ka*  < 

^  +  ka  then,  while  we  think  we  are  stocking  in  such  a  way  as  to 
achieve  a  given  risk  p;  in  fact  the  true  risk  p*  would  be  greater. 
The  resulting  shortage  "cost"  paid  for  such  a  position  could  be 
disastrous.  Since  neither  position  is  particularly  favorable,  and 
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the  true  parameter  values  are  unknown,  it  is  clear  that  the  most 
precise  estimates  of  ^  and  a  are  desirable;  precise  in  the  sense 
of  minimum  fluctuation  about  the  true  values. 

Trying  to  determine  the  true  risk  incurred  when  exponentially 

smoothed  estimates  |i  and  o  of  ^  and  a,  respectively.,  are  used 

theoretically  requires  the  joint  probability  distribution  of  X,  ^ 

~  ~ 

0.  If  finding  the  distribution  of  o  alone  seems  difficult,  the 
task  of  finding  such  a  Joint  distribution  looms  formidable  to  say  the 
least.  Fortunately,  in  a  simulation  approach  the  parameters,  and 
hence  the  true  theoretical  reorder  level  for  a  given  risk  p,  is 
known.  For  any  procedure  used  to  set  actual  reorder  levels  it  is  then 
possible  to  observe  the  behavior  of  repeated  applications  of  such  a 
procedure.  There  are  many  ways  this  might  be  done.  One  approach  is 
to  observe  the  number  of  times  a  reorder  level  falls  within  so  many 
units  of  the  theoretical  level.  Or,  for  a  given  percentage  P  we 
might  ask  within  how  many  units  of  the  theoretical  reorder  levels 
will  P  percent  of  the  actual  reorder  levels  be  found?  We  have 
done  a  little  of  both  and  summarized  the  findings  for  various  cases 
in  Table  3. 

The  table  is  constructed  is  follows.  First  a  triple  of  parameters 
(pP  P)  is  chosen  where  p  is  the  desired  risk.  From  this  k  io 
determined  from  tables  and  a  theoretical  reorder  level  computed.  For 
example,  if  |i  B  400  0  m  40  and  p  *  .01  thea  from  .01  ■  P(3^400+40k) 

we  determine  k  *  2.326  and  theoretical  reorder  level  is  49.3 . 04.  Then, 
for  each  parameter  triple  Qj  is  the  percentage  (rounded)  of  200 
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reorder  levels  found  to  be  within  1  unit  of  the  theoretical  reorder 
level  first  for  the  level  based  on  exponential  smoothing  fi  +  ko 

^  A  A 

and ,  second;,  for  a  level  set  using  X  +  kc  where  o  is  the  maximum 
likelihood  estimate  of  o  and  X  is  the  average  demand.  Next,  P^q 
represents  the  number  of  units  about  the  theoretical  reorder  level 
within  which  50X  of  the  computed  reorder  levels  were  found.  Again, 
the  first  column  under  P^q  is  for  smoothed  estimates  and  the  second 
for  maximum  likelihood  estimates.  The  last  column  is  a  similar 


computation  for  95X  of  the  reorder  levels. 


Parameter 

Triples 

_ Ql 

p50 

p95 

(100,25,. 01) 

4 

32 

15 

2 

35 

13 

(100, 25,. 11) 

6 

48 

11 

1 

24 

6 

(100,25,. 50) 

8 

19 

6 

2 

17 

5 

(400,40, .01) 

2 

22 

25 

2 

57 

2! 

(400,40,  .11) 

5 

38 

18 

L 

38 

9 

(400, 40,.  50) 

5 

10 

10 

.3 

27 

7 

(600, 60,. 01) 

1 

15 

37 

4 

85 

31 

(600,60, .11) 

5 

25 

26 

.3 

58 

14 

(600, 60,. 50) 

3 

6 

15 

4 

40 

12 

(800,80, .01) 

2 

l5 

35 

5 

103 

46 

(800, 80  ,.11) 

2 

If 

28 

3 

69 

30 

(800, 80,. 50) 

2 

2.3 

21 

2 

56 

23 

(1000, 100,.  01) 

2 

13 

44 

6 

129 

58 

(1000, 100,.  11) 

1 

12 

35 

4 

85 

37 

(1000, 100 ,.50)  ’ 

2 

18 

_ 

2" 

» 

70 

30 

Table  3.  Percentiles  for  Reorder  Levels 
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Many  more  cases  than  those  reported  in  the  table  were,  examined 
but  the  results  are  omitted  for  the  sake  of  brevity.  The  results  in 
these  other  cases  displayed  precisely  the  same  pattern  however. 

Indeed,  the  consistency  of  the  various  cases  is  somewhat  startling. 
Once  again  there  is  not  one  single  choice  of  a  mean,  a  standard 
deviation  or  a  risk  in  which  the  smoothing  procedure  does  not  fall 
significantly  short  of  maximum  likelihood  techniques.  This  was  to 
be  expected  of  course  since  the  variance  of  the  estimators  involved 
determine  the  amount  of  flucutation  from  sample  to  sample. 

One  other  interesting  observation  in  Table  3  is  the  monotonic 
nature  of  P^q  and  P^.  Thus,  while  it  always  takes  more  units  of 
distance  (roughly  twice  as  many)  to  include  93X  of  the  values  using 
smoothing  versus  maximum  likelihood,,  the  comparison  at  the  50X  level 
is  even  more  striking.  In  addition  to  requiring  roughly  ten  times  as 
many  units  to  cover  50?  of  the  observations  for  smoothing,  the 
number  of  units  for  tne  same  percentage  using  maximum  likelihood  is 
comfortably  small. 

Another  way  to  view  the  results  is  to  try  and  estimate  the  true 
risk  incurred  when  a  given  procedure  is  followed.  One  way  to  do  this 
is  to  imagine  that,  after  200  observations,  the  estimates  say  l* 
and  coincide  with  the  parameters  and  then  compute  P(]^*+k<7*) 
as  e  measure  of  the  true  risk.  With  a  sample  size  of  200,  the 
asymptotic  result  of  supposing  that  estimates  and  parameters  are  the 
same  should  be  fairly  sound.  Of  course  the  point  is  that  how  sound 
such  an  estimate  is  depends  upon  the  precision  of  the  estimate  once 
again.  In  any  event  in  practice  this  is  essentially  what  is  done  and 
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200  is  a  figure  quite  a  bit  larger  than  'he  typical  number  of  periods 
for  which  historical  data  are  aaintained  for  making  such  estimates 
within  NavSup. 

In  Table  4  our  results  are  summarized  as  follows.  Again,,  certain 
triples  of  parameters  p,  cr  and  p  are  selected.  This  time,  for 
each  pair  (p,  a)  down  a  column  five  different  risks  are  associated 
across  a  row.  For  each  such  combination  an  actual  reorder  level  is 
determined  in  two  ways.  First  of  all  k  is  determined  and  then 


reorder  levels  fixed  at  p  +  ka  and  X  +  k o .  After  the  200th 
observation  p*  ■  P(X>p*+kC*)  is  computed  and  entered  in  each  cell 

it 

(suitably  rounded),  first,  for  (p  ,o  >  «  (p„a)  and  second  for 
(p*,ff*)  ■  (X.cr)  For  lack  ot  a  better  name*  p,;  is  called  the  actual 
risk.  In  this  way,  the  actual  risk  incurred  using  smoothing  may  be 
compared  with  maximum  likelihood  procedures. 


(p^y^ 

.01 

.05 

.11 

.25 

.50 

( 100  j,  25) 

.01 

.01 

.05 

.05 

.1.3 

.11 

.30 

.25 

.60 

.52 

(400,40) 

.01 

.01 

.05 

.05 

.15 

.11 

.30 

.25 

.60 

.52 

(500,50) 

.01 

.01 

.05 

.05 

.  13 

.11 

.30 

.25 

.60 

.52  ! 

(600.60) 

.01 

.01 

.05 

.05 

.13 

.11 

,30 

.25 

.60 

.52 

(800,80) 

.06 

.01 

.18 

.05 

.31 

.12 

.52 

.26 

.76 

.52 

(900,90) 

.12 

.02 

.20 

.0" 

.27 

.13 

.37 

.27 

.50 

.51 

(1000,100) 

.06 

.01 

.18 

.05 

.31 

.12 

.52 

.26 

_ 

.77 

.52 

Table  4,  Acn.al  Risks 

The  results  are  fairly  self-explanatory.  The  results  are  even  fairly 
consistent  in  the  sense  that  for  most  parameter  choices  the  asymptotic 
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results  are  quite  the  same.  For  large  means  and  variances  more 
variation  begins  to  show  up.  We  see  that  for  small  values  of  p,  both 
procedures  are  roughly  equivalent.  That  is  to  say,  after  200 
predictions,  the  actual  risk  1s  fairly  close  to  the  true  risk.  But 
for  large  or  even  moderate  values  of  p  it  appears  that  smoothing 
consistently  produces  a  larger  actual  risk  than  maximum  likelihood. 
Both  are  larger  than  the  true  risk  but  the  maximum  likelihood 
procedure  is  much  closer  to  the  true  value  and  is  far  more  consistent. 


Risk 

Ji.' 

M*  +  ko 

P  +  ko 

x  +  ka 

.01 

51 

38 

44 

.02 

52 

36 

43 

.03 

52 

40 

41 

.04 

53 

39 

40 

.05 

54 

39 

39 

.06 

57 

38 

37 

.07 

5? 

40 

37 

.08 

52 

40 

34 

.09 

^2 

40 

35 

.10 

51 

4] 

34 

Table  5.  Percentage  Overestimates 
Finally,  supposing  the*-  it  is  far  more  serious  to  overestimate 

the  reorder  level  (resulting  in  excess  slock  for  a  given  risk)  than 
to  underestimate  it.,  the  data  were,  analyzed  to  see  the  percentage  of 
time  the  reorder  level  was  overestimated  in  200  trials.  For  this 
illustration  only  disti  ib-<*  i  on  was  used  namely  normal  with 
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mean  50  and  standard  deviation  10.  Again,  for  various  risks 
(this  time  no  more  than  10%)  reorder  levels  are.  set  using  +  k0 

D  « 

and  x  +  ka.  This  time  we  examined  a  thtTd  alternative  a  mixture 

,U. 

of  the  two  procedures;,  given  by  jji  +  k<J.  The  results  are  summarized 
in  Table  5  and  the  entries  are  the  percentage  of  times  that  the  given 
method  overestimates  the  theoretical  reorder  level. 

It  is  by  this  Ime  perhaps  not.  surprising  to  note  that  smoothing 
consistently  overestimates  the  reorder  level.  Even  the  mixed  procedure, 
whereby  smoothing  is  used  to  estimate  the  mean  and  a  is  estimated  by 
maximum  likelihood,  is  an  improvement  over  a  MAD  estimate  of  o. 

Compared  with  straight  maximum  likelihood,  it  is  better  for  small 
values  of  p  and  then  begins  to  fall  off  as  p  increases. 

5 •  CONCLUSIONS  AND  RECOMMENDATIONS 

In  drawing  conclusions  it  would  be  well  to  reiterate  the  tentative 
nature  of  simulation  results.  Without  the  exact  probability  distri¬ 
bution  for  the  exponentially  smoothed  estimators  treated  in  this 
report—particularly  MAD— no  definite  statement  can  be  made  about  the 
efficiency  of  this  method  over  other.  Bu'  the  simulation  results 
display  a  certain  consistent v  that  'eitainlv  lends  support  to  the 
claim  that  maximum  likelihood  methods  would  be  superior  tor  those 
purposes  to  which  smoothing  is  presently  employed.  Certainly  the 
evidence  is  sufficient  t"  warrant  t-  rther  examination  and  comparison 
with  alternatives.  In  spite  of  the  difficulties  involved  in  the 
derivation  of  the  probability  distributions  when  exponential  smoothing 


Is  employed  as  a  forecasting  fool,  it  is  our  recommendation  to  continue 
the  theoretical  search  for  those  distributions  or  at  least  the  moments 
that  are  involved.  In  particular  it  would  be  helpful  to  establish 

1*3 

that  the  variance,  of  o  is  bounded  below  by  a  positive  quantity  as 
seems  to  be  indicated  by  the  simulations.  If  this  were  true  then  the 
estimate  will  always  fluctuate  abc.it  the  true  o  regardless  of  how 
many  time  periods  are  considered. 

Another  reason  for  continuing  research  in  this  area  is  the  fact 
that  aven  the  simulation  studies  were  only  carried  out  for.  the  normal 
case  (and  the  isolated  Poisson  cases)  with  constant,  mean  and  variance. 
This  is  a  strong  assumption  and  entails  supposing  that  no  trend  is 
present.  It  is  highly  recommended  that  a  similar  analysis  in  the 
presence  of  trend  be  carried  out.  The  trend  test  and  indicator 
presently  used  by  NavSup  as  outlined  in  [,4]  agai*  employs  exponential 
smoothing  as  a  basic  tool  for  analyzing  the  data.  If  the  lack  of 
efficiency  indicated  in  this  st~dy  carries  over  to  that  case  as  well 
then  the  consequences  should  be  even  score  serious.  Indeed,  there  is 
reason  to  believe  that  with  the  methods  presently  employed,  there 
can  be  a  50%  chance  of  rejecting  an  almost  vertical  trend,  that  is 
a  sudden  jump  in  the  demand  pattern.  These  matters  are  certainly 
topical  candidates  for  further  research. 

In  any  or  all  of  the  above,  recoroiendatioris,  it  is  suggested  that 
real  data  be  used  to  analyze  unfortunately  perhaps  in  T«trospect, 
the  reorder  system  had  methods  su:h  as  raax'tuvjm  likelihood  been  used. 
Large  amounts  of  su  h  data  ate  available  tot  a  wide  variety  of 
inventory  items  and  the  pi  ograms  i<>»  prove m  h  data  have  been 
preserved . 

.  t 


6.  APPENDIX 


This  appendix  is  devoted  to  supplying  the  mathematical  details 
for  the  monotonic  nature  of  £  a?  asserted  in  Section  3.  It  should 
be  recalled  from  that  section  that  ^  ,  expressed  as  a  function  of 
the  Poisson  parameter  X  is  given  by 


A 

o 

In  particular,  when  X  ■  n 
ur  given  by. 


2  A 


X  x1 


W 


[X]! 


-  has  a  relative  minimum  with  a  value 

a 


u 

n 


2  /n 


-n  n_ 

nl 


It  is  conjectured  that  this  sequence  of  values  is  monotone  increasing. 
To  see  this  observe  that 


u 


n 


“n+l 


and 


u  , 

log  _JL.  *  (n+%)  log  -— ■  +  1  *  “(n-M{)  log(l-M  +  1  . 
un+l  n  " 

Recalling  the  Taylor  expansion  of  log  (I+x)  (see  [6]  for  example), 
we  have,  for  n  >  l., 


1 

R 


1 


1 

2n  3n3 


R 


where 


L. 


it 


1 


4n" 
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Then^ 


u-  '  1  1  1,1.1,  2nrl 


log  — —  B  "[l  +  — i 


U 


n+1 


2*  2-  4n‘  3n2  6n3 


R]  +  1 


1  _  1  _  2n+l 

I2n2  6n3  2 


R 


But  R  <  0  £or  this  case  so  that 


2n+2  2n +2  ,  ,  1  1 

—  «*—  |r!  *”+  — 

4nJ  8iT 


Hence, 


log 


n+1 


1  11.1 

.  ‘  71  71  7? 

12n  6n  4n  8n 


( i'n4"  2n-  3) 


n/,_ 

£,40 


Now  for  x  l  2,  the  polynomial  2x4  -  2x  -  1  Is  positive  so  that 


u  u 

log  — JL  <  0  and  hence  — 3-  <  1 
un+l  un+l 


In  other  woids. 


u  <  ^  .  n  « 

n  n+1 


which  is  what  was  to  be  verified. 
From  [7]  (page  h+) t 


1 


e"(n+w/2i  (r.+Vn+4  ,  fMp,V  <  r!  <  /in  {-.*:**  . 


we  have 


Since  f  (n+fc)  ■ 


Thus  the  sequence 
the  claim  on  page 


2 


1 


ft  <  f  (n+fc)  <  ft  e  24(D^ 
TT  “  IT 


of  relative  maxima  of  ^  converges  to  ft  supporting 

A  " 

7  that  *  then  converges  to  this  limit. 
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